Nature Plants

Plant-specific RNA Polymerase V (Pol V) transcribes noncoding RNAs in the RNA-directed DNA methylation pathway, thereby influencing gene expression and genome stability by controlling de novo DNA methylation. However, the mechanisms governing precise chromatin localization and transcriptional activities of Pol V remain elusive. Here we show that Pol V localization is spatially constrained by the chromatin regulators microrchidia (MORC) and MORPHEUS' MOLECULE 1 (MOM1). MORC and MOM1 promote Pol V occupancy at sites near active chromatin, whereas their loss leads to redistribution of Pol V into CMT3-enriched heterochromatin, accompanied by noncoding RNA transcription, small RNA production and DNA methylation. Our findings reveal a combinatorial model in which recruitment, spatial constraint and DNA methylation feedback collectively define Pol V chromatin distribution and epigenetic function.

Oil palm (Elaeis) is the primary source of global vegetable oil. Interspecific hybrids of Elaeis exhibit pronounced heterosis by integrating two distinct subgenomes into a single nucleus, effectively combining the high yield of African oil palm (E. guineensis) with the high unsaturated fatty acid content and disease resistance of American oil palm (E. oleifera). However, the genetic basis underlying heterosis is still unclear. Here, we combine phased genome assembly, comparative genomics, evolutionary genomics and haplotype-aware transcriptomics to unravel the genetic architecture of heterosis of hybrid oil palm. We assemble the highly heterozygous F1 genome ('Reyou 40', 3.75% heterozygosity) into a complete 1.73 Gb T2T haplotype (HapG) and a 1.84 Gb near-T2T haplotype (HapO with17 gaps). Despite 91.56% sequence identity, HapG and HapO diverged in LTR-RT occurrence and PAV affected genes, showing complementary biases in lipid metabolism and stress responses, respectively. Evolutionary genomics revealed that ancient WGDs preserved the palm family. Whereas lineage-specific lipid-related gene expansions in oil palm. Six ancient introgressed regions (~64 Mb) in HapG were reshaped by transposable elements and tandem duplication, showing an enrichment of genes related to resistance and lipid metabolism. Transcriptomically, 82.2% of allelic gene pairs maintained balanced expression, accompanied by parental functional complementarity and dosage buffering, revealing a potential regulatory basis for coordinating parental genetic differences in the hybrid genome. These haplotype-resolved genomic resources offer vital targets for understanding heterosis and accelerating oil palm molecular breeding.

Anterograde and retrograde signalling establish bidirectional communication between the nucleus and chloroplasts. Retrograde signals from chloroplasts regulate nuclear gene expression while anterograde signals from the nucleus coordinate chloroplast development and maintain cellular homeostasis. How this bidirectional signalling framework extends beyond locus-specific regulation to shape the global spatial organization of nuclear chromatin across tissues remains unclear. Although anaesthesia can alter chromatin organisation, the role of chloroplast dysfunction in these changes remains unclear. Here, we investigate how chloroplast dysfunction and anaesthesia influence euchromatin and heterochromatin organisation in Solanum lycopersicum seedlings across tissues with contrasting photosynthetic competence. Using confocal and super-resolution radial fluctuation (SRRF) imaging with quantitative multiparameter analysis, we identify distinct, tissue-specific chromatin responses to chloroplast disruption and anaesthesia. Notably, anaesthesia induces distinct spatial chromatin changes across tissues that are independent of chloroplast dysfunction, suggesting a direct nuclear response to anaesthesia rather than a chloroplast-mediated retrograde effect. These findings highlight chromatin topology as a potential quantitative biomarker of cellular disruption and provide a framework for investigating anterograde chloroplast-nucleus coordination and stress-responsive nuclear organisation in plants.

Antibiotic-induced disruption of plant-associated microbiomes has the potential to alter host health beyond the directly exposed generation, yet whether the effects of dysbiosis are transmitted through the seed microbiome remains unknown. Here, we investigated the intergenerational impacts of streptomycin-induced dysbiosis in tomato (Solanum lycopersicum), demonstrated that seed microbiome transfer (SMT) restores progeny microbiome function and disease resistance, and characterized the underlying physiological and genetic mechanisms. Parental streptomycin exposure altered the composition of progeny rhizosphere bacterial communities, reduced expression of defense-associated genes, and increased susceptibility to Xanthomonas perforans. Suppression of immune gene expression was strongly associated with increased disease severity, indicating that parental dysbiosis impaired progeny plants ability to mount effective immune responses. Transfer of the seed microbiome from healthy plant donors partially restored rhizosphere community composition, reduced disease severity and recovered defense gene expression of three genes. Together, our findings demonstrated that antibiotic exposure microbiome disturbance generates intergenerational legacy effects that influence plant immunity and disease susceptibility and seed microbiome transfer can counteract this dysbiosis across generations.

Zinc is an essential micronutrient at low concentrations, yet it becomes toxic at slightly higher ones. This is exploited by plants as an effective defensive strategy. However, the molecular components that are involved zinc-mediated immunity remain poorly defined. Here, we show that mixed-linked {beta}-1,3/1,4-glucans naturally occurring in microbial and grass cell walls and used as an agrobiological solution, trigger zinc accumulation in the Arabidopsis apoplast and upregulate the expression of the zinc transporters HMA2 and HMA4. This response occurs independently of salicylic acid, jasmonic acid and ethylene-mediated signalling pathways, but it requires the LysM receptor kinases CERK1, LYK4 and LYK5, indicating a specific pattern triggered immunity-associated mechanism. We further demonstrate that hma2hma4 mutants display constitutive activation of a broad set of defence-related genes, yet this transcriptional reprogramming is insufficient to confer resistance against the necrotrophic fungus Plectosphaerella cucumerina BMM. Moreover, metabolomic profiling highlights the contribution of specialized metabolites to this defective defence output. Altogether, our findings reveal that zinc-mediated toxicity constitutes a defence mechanism integrated into the immune response triggered by specific microbial or damage associated molecular patterns.

MicroRNAs (miRNAs) are key post-transcriptional regulators of development, stress response, and secondary cell wall formation in woody plants, yet annotations for Eucalyptus grandis, the world's most widely planted hardwood, remain fragmented across studies using incompatible discovery pipelines and filtering criteria. Here we present the Eucalyptus MicroRNA Archive (EMA), a curated, locus-resolved database integrating three independent small RNA sequencing datasets spanning vegetative tissue, somatic embryogenesis, and mechanically induced tension wood formation. Applying annotation criteria aligned with current plant miRNA standards, EMA catalogs 99 curated miRNAs (31 previously described, 68 novel) organized into 34 family-level groupings under a three-tier confidence system, known-reference-supported, multi-study replicated, or single-study, that preserves study-of-origin and sample-level evidence for every entry. Cross-study comparison showed that only 9 of 99 entries (9.1%) were independently supported by all three datasets, supporting an evidence-tiered rather than binary annotation scheme. Target prediction against the E. grandis transcriptome yielded 1,773 miRNA-target interactions spanning 764 loci, integrated into a combined miRNA-target and protein-protein interaction network. This network resolved into functionally coherent, mutually isolated clusters, including an miR482-associated NBS-LRR/TIR disease-resistance hub with a substantial translational-repression component, alongside modules enriched for ribosome biogenesis and translation, DNA replication, and nitrogen and carbohydrate metabolism. EMA is publicly accessible through an interactive web dashboard, with all curated data, source code, and analysis scripts openly available, providing a reproducible, extensible framework for E. grandis miRNA research and a template for similarly structured resources in other non-model woody species.

Aegilops peregrina is a wild allotetraploid wheat wild relative and an important source of genetic diversity for stress tolerance and agronomic traits. Here, we report a subgenome-resolved, chromosome-scale reference genome assembly of a drought tolerant and stem rust resistant Ae. peregrina accession PI 604178 generated using PacBio HiFi and Hi-C sequencing. The 10.13 Gb assembly contains 98.81% of sequence anchored to 14 pseudomolecules representing the seven S and seven U chromosomes, with contig and scaffold N50 values of 25.84 and 746.48 Mb, respectively. The assembly achieved a consensus quality value of 74.61, 97.83% k-mers completeness, and 99.9% BUSCO completeness. LTR Assembly Index values of 20.43 and 18.79 for the S and U subgenomes, respectively, further supported high continuity across repeat-rich regions. Repetitive elements comprise 85.93% of chromosome-anchored assembly. We annotated 59,910 high-confidence protein-coding genes, with comparable gene representation across the two subgenomes. This reference genome provides a high-quality genomic framework for comparative analyses, characterization of important loci regulating agronomic and resilience related traits, and sequence-guided exploitation of Ae. peregrina allelic diversity for wheat improvement.

Extrafloral nectaries (EFNs) are a widespread plant defense mutualism trait and are highly convergent, appearing in hundreds of plant lineages worldwide. Here we investigate a report of possible EFNs in Erythranthe angulosa, a recently-described California wildflower. We integrate field observations, insect bioassays, an induction experiment, and microscopy to test for signatures of EFN function, finding no evidence that the distinctive axillary swellings produced by E. angulosa function as EFNs. We also uncovered two distinct morphs at the type locality of E. angulosa that diverge in the number of axillary swellings produced, as well as other shoot architecture traits such as stem thickness, leaf size, and branch number. Although the axillary swellings appear to not function as EFNs, they remain a compelling morphological variant within the yellow monkeyflowers that may perform storage or another unknown function.

Blumeria graminis f. sp. avenae (Bga), the causal agent of oat powdery mildew, is one of the most host-specialized members of the B. graminis species complex. Despite its agricultural importance, the lack of a high-quality reference genome has limited studies of host specialization, virulence evolution and comparative genomics in this pathogen. Here, we generated the first chromosome-scale genome assembly of Bga using an integrative approach combining long- and short-read sequencing, Hi-C scaffolding and transcriptome data. The Bga genome exhibits hallmark features of powdery mildew fungi, including extensive repeat content and low gene density. Comparative analyses revealed that genome expansion is primarily associated with historical transposable element proliferation rather than recent transpositional activity. Genome organization is consistent with a functionally stratified "one-speed" model, in which genes associated with pathogenicity, including predicted effectors and infection-responsive genes, are preferentially located in transposable element-rich regions characterized by reduced synteny conservation and extended intergenic spaces. In contrast, conserved genes are concentrated in compact genomic regions and maintain strong syntenic conservation across cereal-infecting formae speciales. Hi-C analyses demonstrated a highly structured chromatin architecture and revealed genome organization patterns associated with infection-related gene expression. Comparative genomic analyses indicated that host specialization in Bga is driven by localized diversification of a relatively small subset of genes rather than large-scale genome restructuring. These results provide the first high-quality genomic resource for Bga and offer new insights into the evolutionary mechanisms underlying host specialization in powdery mildew fungi.

The integration of digital technologies for high-throughput field phenotyping is critical for accelerating crop improvement in agriculture. However, extracting traits from remote sensing data remains constrained by fragmented workflows, manual intervention, and limited interoperability among existing tools, resulting in delays that hinder timely biological insight and decision-making. To address these challenges, we present PhenoStream (Phenotyping Streaming), a scalable, end-to-end cyberinfrastructure designed to automate the full lifecycle of aerial imagery-based phenotyping, from data acquisition to plot- and genotype-level inference. The framework integrates automated data ingestion from distributed field sites, geospatial processing, and AI-enabled trait extraction within a unified, user-accessible graphical interface. Its modular and extensible architecture supports adaptable trait modeling and seamless integration of new data sources, enabling deployment across diverse crops, environments, and experimental designs. We demonstrate the system across a large multi-location field trial network of bioenergy crops, where it enables high-throughput characterization of spatiotemporal growth dynamics, genotype-by-environment (GxE) interactions, and predictive modeling of key agronomic traits. By significantly reducing processing latency and manual effort, the platform facilitates near-real-time analysis and reproducible workflows. This work establishes a generalizable and scalable pathway for operationalizing very-high-spatial resolution aerial phenotyping in agricultural research. By bridging data acquisition and analytics, the end-to-end cyberinfrastructure provides a foundation for integrating heterogeneous and unstructured data streams--including remote sensing, environmental, and management data--toward data-driven decision making in agriculture.

Genes encoding novel protein sequences are a ubiquitous feature of genomes. They fuel molecular and cellular evolutionary innovations and frequently contribute to species-specific characteristics. We are now unravelling the processes by which they originate, including de novo from noncoding sequences and through extreme divergence, yet how much and what types of novel proteins evolve through each process is still unclear Does the mechanism of origination shape the structural and functional potential of the resulting proteins? Here, we conducted a broad computational investigation of genetic and protein novelty at the scale of the entire subphylum of Saccharomycotina yeasts. We detected more than 5,000 robust de novo genes across 332 species and compared them to more than 10,000 novel genes resulting from extreme sequence divergence, revealing two distinct modes of evolution of novelty. A remarkable 40% of de novo proteins are predicted to localize to mitochondria compared to only 15% of divergent, with the latter also being substantially longer and more disordered. A detailed analysis of conservatively predicted tertiary structures of novel proteins shows that "invention" of novel folds can happen through both processes but is more likely to occur de novo. We also illustrate cases of evolutionary "re-invention" of existing protein folds from non-coding sequences. Our work deepens our understanding of the origins and importance of novel proteins opening new directions for further structural and functional characterization.

Proliferating cells maintain their cytoplasmic density within a narrow range but deviate when entering quiescence or experiencing stress, suggesting active regulation. The mechanisms driving these density adjustments and their impact on cellular function remain unclear. Here, we demonstrate that the conserved cAMP-activated protein kinase A (PKA) is a key regulator of cytoplasmic properties. Inactivation of PKA leads to a drastic increase in cytoplasmic dry mass density and reduced diffusion that depends on the environmental stress response (ESR) transcription factors Msn2/4. This change is mediated by the accumulation of glycogen and trehalose, which have opposing effects on intracellular diffusion. Importantly, the accumulation of these carbohydrates confers stress resistance in distinct ways and independently of their roles as energy sources. Our findings highlight the importance of the biophysical properties of the cytoplasm in stress resistance and the role of glycogen and trehalose in regulating these properties.

Plant cold acclimation emerges from coordinated adjustments in photosynthesis, primary metabolism, and intracellular carbon allocation. Yet, the regulatory role of subcellular metabolite compartmentation in natural variation of cold acclimation remains insufficiently understood. Here, we investigated four Arabidopsis thaliana accessions grown either individually or in bulk to determine how growth configuration and genotype shape the metabolism of sugars and organic acids during cold exposure. Using non-aqueous fractionation, we quantified plastidial, cytosolic, and vacuolar sugar pools alongside whole-cell carbohydrates, organic acids, enzyme activities, photosynthetic parameters, and stress markers. A neural-network classifier revealed that subcellular sugar distribution together with sugar amounts and organic acids provided the strongest discriminatory power among accessions, surpassing photosynthetic traits and enzyme activities. Our findings demonstrate that natural variation in cold acclimation is strongly determined by genotype-specific subcellular metabolite architectures, and that the cultivation strategy modulates these intracellular signatures. We conclude that subcellular compartmentation of metabolites represents a cellular control layer for natural variation of cold acclimation and resilience in Arabidopsis thaliana.

Circadian rhythms in transcription are facilitated by well-defined genetic circuits, but how molecular clocks drive daily rhythms in mammalian physiology is poorly understood. The mechanistic target-of-rapamycin (mTOR) complex integrates daily systemic and circadian intracellular timing cues for input into the cellular timekeeping machinery. Here we demonstrate that mTOR is a major clock output pathway whose activity is required for most daily variation in cellular and organismal physiology, with PERIOD2 shown to interact directly with mTORC1. Acute mTOR inhibition abolishes functional rhythms in cells and most daily variation in mouse liver physiology. mTOR activity is not required for clock protein or locomotor rhythms, indicating that mTOR is not part of the cellular or central circadian timekeeping mechanism. In the forebrain, mTOR activity is required for most detectable daily rhythms in protein abundance and phosphorylation; however, the daily architecture of the sleep/wake cycle is remarkably preserved in mice and zebrafish under mTOR blockade, with a significant increase in wakefulness. Clock outputs in Arabidopsis (plant) and Neurospora (fungus) are also more sensitive to mTOR inhibition than core clock mechanisms indicating evolutionary conservation of mTOR as a circadian effector. We conclude that most but not all daily physiological rhythms in mammalian cells and tissues depend on rhythmic regulation by the mTOR pathway.

Life began with the emergence of a molecule that could replicate its own genetic material, a task plausibly mediated by an RNA-dependent RNA polymerase ribozyme. Here, we present the structure of such a polymerase ribozyme, bound to RNA substrates comprising the template, primer, and nucleoside triphosphate (NTP) analog. The structure reveals how directed evolution shaped flanking elements around a highly conserved catalytic core derived from the ancestral class I ligase ribozyme. Each element serves as a functional module, positioning the primer-template duplex and incoming NTP within the active site of the enzyme. This emergent domain organization is remarkably similar to the "right hand" configuration of polymerase proteins, suggesting a common functional form for copying nucleic acids, regardless of biopolymer catalyst.

Protein synthesis is dynamically regulated to control cell growth, differentiation, and stress responses. Recent single-cell sequencing methods can map ribosome positions on individual transcripts, but cannot capture the global translational states that coordinate protein synthesis across the transcriptome. In contrast, methods that measure the global translational landscape, such as polysome profiling and cryogenic electron tomography, lack either single-cell resolution or throughput. Here we introduce SCISSOR (Single-Cell Inference of Structural States of Ribosomes), a strategy that infers global translation activity in individual cells from the differential protection of ribosomal RNA (rRNA) against nuclease digestion. By integrating these protection signatures with the structure of the ribosome, SCISSOR resolves multiple ribosomal states and quantifies their abundance across thousands of individual cells. Applying SCISSOR reveals systematic variation in global translation across the cell cycle in human cells, as well as during the differentiation of murine intestinal stem cells into distinct epithelial lineages. These findings uncover principles of global translational regulation that are invisible to transcriptomic or ribosome-profiling assays, establishing a framework for studying global translation control at single-cell resolution.

All living organisms rely on the movement of ions across cell membranes as the fundamental physical basis of their internal energy and signaling, and plants are no exception. Plants perceive, integrate, and respond to environmental stimuli through electrical signals, classified as action, variation, and system potentials, that are coupled with calcium waves, reactive oxygen species, and hydraulic and hormonal changes to coordinate whole-organism responses despite the absence of a nervous system. Yet most studies characterize these signals using a single feature, such as amplitude or spike duration, in a single tissue, an approach that cannot establish how such signals correspond to the underlying ionic activity, mobility, and structural complexity of the signaling environment, or how this correspondence varies across organs. Here, we correlate plant bioelectrical signals with potential ionic energy flow using a multi-domain framework, combining discrete spike events, continuous waveform properties, spectral composition, and signal complexity applied to leaf, stem, and root recordings from tomato (Solanum lycopersicum) exposed to different stimulus. Electrical activity with increased stimulus strength, likely reflecting increased ionic flow, with the root showing the largest response. This suggests plant electrical signaling works as a distributed, ion-based information system, useful for stress monitoring and bio-inspired sensor design.

Accelerating photoprotective regulation to improve carbon assimilation is a promising strategy to increase crop productivity. Although rapid non-photochemical quenching (NPQ) relaxation has been validated as a target through metabolic engineering, it remains unclear whether conventional breeding has improved this trait. Here, we investigated whether more than a century of soybean breeding enhanced NPQ relaxation alongside light-saturated carbon assimilation and seed traits. We evaluated a historical panel of 24 soybean genotypes across vegetative and reproductive developmental stages by integrating NPQ relaxation, gas exchange parameters, xanthophyll-cycle pigment profiles, expression of key photoprotective genes (VDE, PsbS, and ZEP), seed number and seed weight. NPQ relaxation parameters were not consistently associated with genotype release year, seed number, or seed weight at either developmental stage. The only exception was the amplitude of the rapidly relaxing NPQ component (AqE), which was negatively correlated with all three variables during the reproductive stage. In contrast, genotype release year was positively associated with maximum net CO2 assimilation rate (Amax), maximum carboxylation rate of Rubisco (Vcmax), maximum electron transport rate (Jmax), seed number, and seed weight, while Amax and Vcmax were positively correlated with seed number and seed weight. These findings indicate that the greater photosynthetic capacity of modern genotypes was not accompanied by faster photoprotective response. Thus, photoprotective regulation has not kept pace with gains in photosynthetic capacity under field conditions. We conclude that rapid NPQ relaxation remains an important target for synchronizing photoprotection with the high photosynthetic capacity of modern soybean lines.

Drug resistance emerges repeatedly in outbreaks of Candida fungal pathogens, but little is known about its origins or persistence. Here, we investigated the evolutionary processes shaping echinocandin resistance in Candida auris, a globally emerging and predominantly clonal fungal pathogen. Genome-wide association across over 600 isolates identified mutations in the {beta}-1,3-glucan synthase gene FKS1 as the most significant driver of resistance to an echinocandin drug. Ancestral reconstruction of this population traced shared resistance mutations among small groups typically consisting of 2-3 closely related isolates, but clusters could include up to 16 isolates. Nearly all resistant clusters consisted of isolates collected in the same year and region, consistent with local transmission. To further examine population-level selection, we measured adaptive signatures in FKS1 and the highly diverged paralog FKS2 across 22,000 genomes. This revealed excess nonsynonymous polymorphisms in FKS1, primarily due to independent, recurrent mutations at resistance hotspots, consistent with parallel evolution and incomplete fixation of adaptive alleles. In FKS2, there is no evidence of hotspots and little support for diversifying selection. Together, these results indicate that resistance mutations emerge under strong genetic constraint, with adaptation restricted to only one FKS homolog and predominantly at mutational hotspots.

Protein language models organize sequence and structure at scale, but a global representation of how proteins respond to mutation remains lacking. We present RegimeFormer, a large protein perturbation model coupled to RegimeAtlas, constructed by harmonizing and indexing 202,556,313 non-redundant protein sequences across the tree of life. A diversity-preserving one-million-protein subset provides the high-resolution training and inference layer, with 995,995 proteins yielding residue-level summaries across 407,048,356 residues and substitution-specific predictions available on demand. Across experimental deep mutational scanning, molecular benchmarks, structural confidence and evolutionary constraint, RegimeFormer identifies reproducible protein-level perturbation regimes that organize residue fragility, adaptability and predictive uncertainty. Regime conditioning improves substitution-specific prediction, with the largest relative gains under unseen-protein, unseen-family and low-homology evaluation. RegimeFormer-derived molecular priors further improve downstream transcriptomic and drug-response modelling. Together, RegimeFormer and RegimeAtlas provide a scalable framework for mapping, predicting and querying protein perturbation landscapes across global sequence space.